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The many-body localization transition is a dynamical quantum phase transition between a local¬ 
ized and an extended phase. We study this transition in the XXZ model with disordered magnetic 
field and focus on the time evolution following a global quantum quench. While the dynamics of the 
bipartite entanglement and spin fluctuations are already known to provide insights into the nature 
of the many-body localized phases, we discuss the relevance of these quantities in the context of 
the localization transition. In particular, we observe that near the transition the long time limits of 
both quantities show behavior similar to divergent thermodynamic fluctuations. 
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I. INTRODUCTION 


Many-body localization (MBL) occurs when Anderson 
localizatiorP persists in the presence of interactions. In 
the pioneering work of Basko, Aleiner and Altshuler,!^ the 
localized phase was shown to be perturbatively stable to 
small interactions. This work quickly opened up a new 
field and many intriguing properties of this new phase 
were explored: (i) due to the lack of transport, MBL 
systems do not thermalize,!^ (ii) at finite energy densities 
the localization of domain walls allows stabilizing quan¬ 
tum and topological order which would otherwise rnelt,!^ 
and (iii) following a global quantum quench, MBL phases 
have a characteristic logarithmic growth of entanglement 
as a function of time.^^J®^ On the experimental side, first 
progress has been made in realizing such systems: In 
Ref. [11] the effect of localization was observed in a cold 
atom experiment where a charge density wave failed to 
relax in the localized phase. By measuring I — V char¬ 
acteristics of amorphous iridium-oxide. Ref. [12] has pro¬ 
vided evidence for a finite temperature insulator where 
the MBL mechanism might be at play. 

As the MBL transition occurs in eigenstates at finite 
energy densities instead of just the ground state, this 
transition is a dynamical quantum phase transition.^ 
Many aspects of this transition from an MBL phase to an 
extended one are still not fully understood. An interest¬ 
ing feature of the transition is that, in principle, the crit¬ 
ical disorder strength depends on energ y density, yield¬ 
ing a so-called many-body mobility edgeW^^^^^^ Novel 
real space renormalization group methods have been de¬ 
veloped in which this tran sition is given by an infinite 
randomness RG fixed point 

In this work, we consider the anti-ferromagnetic spin- 
1/2 XXZ chain and study the time evolution of the en¬ 
tanglement as well as the bipartite fluctuations follow¬ 
ing a global quench. We focus on the evolution of the 
probability distribution and show that the standard de¬ 
viation of the long-time limit can be used to detect the 
MBL transition. The observations made for the bipartite 
fluctuations are particularly useful for an experimental 
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FIG. 1. A schematic representation of a ID spin chain sub¬ 
divided into two parts A and B, which are used to calculate 
the entanglement entropy and bipartite fluctuations. 


detection of the transition in cold atomic systems. We 
furthermore discuss the behavior after a very long time 
following the global quench in comparison to that in the 
thermal state and in the diagonal ensemble. 

This paper is organized as follows: In section |n| we de¬ 
scribe the model and briefly mention some of its proper¬ 
ties. In section [in] we describe the global quench protocol 
followed by the description of the behavior of entangle¬ 
ment and bipartite fluctuations. We then compare the 
long time behavior to that in the thermal state and di¬ 
agonal ensemble. We present the quench results for the 
non-interacting case next and finally conclude by provid¬ 
ing a summary and outlook in section [LV] 


II. MODEL 

We consider the anti-ferromagnetic spin-1/2 XXZ 
model on a one-dimensional chain in the presence of a 
disordered 2 :-directed magnetic field. The Hamiltonian is 
given by 

L-l L 

H = J + + + X hiSt, (1) 

i=l i=l 

where J > 0 is the anti-ferromagnetic coupling strength 
between neighboring spins, A is the anisotropy parame¬ 
ter and hiA are the uncorrelated random external fields. 
Throughout this paper we consider the case J = A = 
1 (except for the non-interacting case when A = 0) 
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FIG. 2 . (color online), (a-b) Time evolution of the entangle¬ 
ment entropy averaged over disorder realizations, for different 
disorder strengths 77 = 1,10 and system-sizes L. The mean 
entanglement saturates after an initial growth for both weak 
(a) and strong disorder (b). In the case of strong disorder it 
shows a logarithmic growth over several decades before satu¬ 
rating to a much lower value compared to the weak disorder 
case. To highlight the qualitative difference between strong 
and weak disorder we also show the evolution of the distribu¬ 
tion of entanglement, color scale, for L = 12 in (c, d). 


and choose hi from a uniform distribution [— 77 , 77 ] . 11 ^ 
The Hamiltonian 0 provides a simple mod el to study 
the MBL phenomena numericallyl ^ l ^^^ l ^^ l ^^ * This model 
shows a localization transition at 77 c ~ 3.6 at infinite 
temperature cor respo nding to eigenstates in the middle 
of the spectrum.E^Ell it has been argued that MBL sys¬ 
tems have a many-body mobility edge,l^ which was first 
observed numerically in transverse field Ising chain.l^ The 
mobility edge has also been obtained for the XXZ chain 
in Ref. |T4| and spinless fermions in Ref. [151 


III. QUENCH DYNAMICS 

Following Refs. [ 6 ] and [7| we consider a global quench 
starting from a simple product state. In particular, we 
choose the Neel state (a product state of alternating up 
and down spins) as the initial state and study the time 
evolution of the system using exact diagonalization. This 
simulation corresponds to a global, sudden quench in 
which we start from the ground state of Hamiltonian 0 
with an infinite staggered field which is then turned off 
at t = 0. In the following, we perform a detailed anal¬ 
ysis of the time evolution of entanglement and bipartite 
fluctuations. 


FIG. 3. (color online), (a) The change in the distribution 
of saturation entanglement Soo (too = 10 ^®) as a function of 
disorder strength 77 . The mean (b) and standard deviation 
(c) of the saturation entanglement entropy as a function of 
disorder strength for different system sizes. The standard 
deviation behaves like divergent thermodynamic fluctuations 
showing a peak which becomes higher with system size. 


A. Entanglement entropy 

We start by considering the evolution of entanglement 
between two equal partitions of the chain (Fig. [^. For 
pure states, the entanglement entropy is given by the 
von Neumann entropy of the reduced density matrix p 
corresponding to either subsystem. The reduced density 
matrix of the left half of the chain {A) for the state It/;), is 
= Tr^ (|t/;)(t/;|), where we have traced over the degrees 
of freedom of the right-half of the chain (H). The von 
Neumann entropy of the state p^ is then given by 

S —Tr^p^ Inp^, 

where the trace is now over the remaining degrees of free¬ 
dom (A). 

Figure (a-b) show the time evolution of S aver¬ 
aged over disorder realizations for different system sizes 
{L = 6,•••,16) and two different disorder strengths 
(77 = 1,10). We have used 10^ disorder realizations for 
L < 12, 10^ for L = 14 and 500 for L = 16. For weak 
disorder (77 = 1 ) the entanglement shows a fast linear 
growth which then rapidly saturates to a value S'oo- This 
linear growth is due to the spreading of correlations at 
a finite spe ed be fore saturating because of finite size of 
the system. 1^2111 xhe saturation value follows a volume 
law Soo = with a being close to its maximum value 
of o^max = ln( 2)/2 (for a partition of the system into two 
half chains). For strong disorder (77 = 10), the system 
shows a rapid linear growth only for a short duration.!^ 
Then the localization causes the linear growth to termi¬ 
nate and is followed by a slow logarithmic growth for 
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a long time before it eventually saturates for finite sys¬ 
tems. The saturation value Soo still shows a volume law 
but the coefficient a is much smaller than it is in the 
weak disorder case. We note that the duration of log¬ 
arithmic growth increases with increasing disorder and 
system size, and we have chosen large enough time in¬ 
terval to study the saturation properties for our system 
sizes. This logarithmic growth in the localized phase has 
been studied in detail in recent work^^K^l ^^d has been 
explained via an interaction induced dephasing mecha- 
nism.l^ 

To gain further insight into the details of entanglement 
dynamics, we plot the distribution of S in Fig.[^(c-d) at 
different disorder strengths for L = 12, and find that it 
differs strongly between the two cases. In the case of 
weak disorder there is a single peak which broadens and 
shifts to higher value of S before saturating. For the 
strongly disordered case, in contrast, the starting distri¬ 
bution with a single peak splits into a bimodal distri¬ 
bution with two peaks at intermediate times; the larger 
peak being near zero while the smaller and much sharper 
peak is at ln(2). This value of the second peak corre¬ 
sponds to cutting one singl et in the partition between left 
and right half of the chain.^^^^^ The first peak at smaller 
value of S slowly broadens and shifts to slightly larger 
values before its tail merges with the second smaller peak 
and the asymptotic distribution thus has a single peak. 
During this broadening the second peaks stays at ln(2). 
The time interval over which this broadening of the main 
peak happens corresponds exactly to the duration of log¬ 
arithmic growth. 

Figure |^a) shows the asymptotic distribution of en¬ 
tanglement as a function of disorder strength. For very 
weak disorder the entropy distribution is centered around 
a relatively high value. The distribution first broadens 
with disorder followed by the main peak shifting to lower 
values. Both for very weak and for strong disorder, the 
distribution is very narrow. To make these observations 
more quantitative, we also show the mean Soo and stan¬ 
dard deviation SSoc of Soo in Fig- (b-c). The mean 
changes slowly for weak disorder but suddenly decreases 
to very small values as the disorder is increased beyond 
a critical value. The standard deviation behaves similar 
to fluctuations in thermodynamic phase transitions and 
shows a peak which gets higher with system size. The di¬ 
vergence of the standard deviation is due to the fact that 
near the transition small changes in the energy densi¬ 
ties and disorder realizations decide whether the system 
is localized or extended. Thus this quantity is a good 
observable to pinpoint the transition. 


B. Bipartite fluctuations 


We now consider the dynamics of bipartite fluctua- 
tionJ^ following the global quench. While the total 
magnetization of Hamiltonian Q is conserved, the 







FIG. 4. (color online), (a-b) Time evolution of bipartite fluc¬ 
tuations averaged over disorder realizations, for different dis¬ 
order strengths r] = 1,10 and system-sizes L. Just like the 
entanglement, the mean of bipartite fluctuations saturates af¬ 
ter an initial growth for both weak (a) and strong disorder 
(b). However unlike the entanglement entropy it does not 
show a logarithmic growth for strong disorder. We show the 
evolution of the distribution of bipartite fluctuations, color 
scale, for L = 12 in (c, d). 


magnetization of the half-chain 

L/2 

i=l 

fluctuates. We define the bipartite fluctuations as the 
quantum fluctuations of S[^ 2 ^ 

^ = (V’i(5i/2)V) - (V’|5i/2lV’)^- 

Entanglement is difficult to measure experimentally 
though there have been recen t suggestions to observe its 
effects in MBL systemsother hand T can 
be accessed in experiments as follows. There is an exact 
mapping between Hamiltonian 0 and hardcore bosons 
on a ID lattice, and the bipartite fluctuations defined 
above are equivalent to particle number fluctuations in 
the bosonic syst em, w hich can be measured using single 

atom microscopy. F^F9 | 

Figure]^ (a-b) shows the time evolution of the disorder 
averaged T. The mean T grows rapidly and saturates 
at very short time scale both in the extended and local¬ 
ized phases. Unlike entanglement there is no logarithmic 
growth for strong disorder and the bipartite fluctuations 
saturate to a much smaller value almost independent of 
system size.l^ 

We present the distribution of as a function of time 
for weak and strong disorder in Fig. (c-d). The be¬ 
havior at weak disorder strength for both the short and 
long time limit is very similar to that of entanglement. 
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FIG. 5. (color online), (a) The distribution of bipartite fluc¬ 
tuations at long times (too = 10 ^®) as a function of disorder 
strength clearly shows the appearance of second sharp peak 
at 1/4 near the expected transition. The weight of the second 
peak gradually shifts to the main peak with increasing disor¬ 
der, resulting in a decrease of the mean. We show the mean 
(b) and standard deviation (c) of bipartite fluctuations with 
disorder strength for different system sizes. 


In particular, we obtain a peak that broadens as a func¬ 
tion of time until saturation. However we can clearly see 
the difference for strong disorder. The short time evo¬ 
lution is qualitatively similar to entanglement, in that 
T also shows a bimodal distribution with a second peak 
at 1/4.1^ The long time behavior on the other hand is 
very different from that of entanglement. For T the first 
bigger peak does not broaden with time and the second 
peak persists even after a long time. This also corre¬ 
sponds to the absence of logarithmic growth and the dis¬ 
tribution saturates much more rapidly than that of en¬ 
tanglement. The absence of logarithmic growth implies 
that though the many-body wavefunction continues to 
evolve, bipartite fluctuations are not affected by the de¬ 
phasing mechanism and attain their asymptotic values on 
a much smaller time scale. It also implies that though T 
can distinguish between localized and extended phases, 
it is insensitive to the effects of interactions and hence 
can not distinguish MBL from Anderson localization. 

The change in the saturation properties of the bipartite 
fluctuations {Too) with disorder strength also captures 
the transition quite well qualitatively (Fig. |^. A second 
peak at the value 1/4 appears in the asymptotic distribu¬ 
tion near the transition, see Fig.|^(a). We show the mean 
Too and standard deviation 5Too of Too as a function of 
disorder strength 77 in Fig. (b-c). In the MBL phase, 
the individual particles can move around only a short 
distance within some localization length. Therefore the 
particle number fluctuations get contribution only from 
particles which are near the partition and as a result 
should be independent of system size. We find that this 


FIG. 6 . (color online). Asymptotic (too = 10^®) entangle¬ 
ment and bipartite fluctuations compared with their values 
in the thermal state for weak (77 = 1 in a,c) and strong 
disorder (77 = 10 in b,d). We use an unequal partition 
{La — ^,Lb — 8 ) for this comparison as the thermal state 
state predictions are expected to be better if the bath (here 
subsystem B) is larger than the system of interest (A). Long 
time after the quench, both quantities agree with the value 
predicted from the thermal state for L = 12 in the case of 
weak disorder (a,c) while the agreement is very poor at strong 
disorder (b,d). We plot the 2D histograms in the insets which 
show the correlation between the asymptotic and thermal val¬ 
ues in each case. 

is indeed the case for strong disorder strengths. This can 
also be seen in Fig.[^(b) where the time-evolution of T is 
almost independent of system size. Just as in the case of 
entanglement the standard deviation behaves like ther¬ 
modynamic fluctuations with the peak becoming more 
pronounced with system size. 

C. Comparison to thermal state 

The notion of thermalization in a closed quantum sys¬ 
tem implies that a generic system would eventually re¬ 
lax and its asymptotic behavior can be described by a 
thermal state with the same energy density. In the case 
of integrable systems, one needs to take into account all 
conserved quantities instead of just the energy,!^ however 
Hamiltonian 0 is not integrable for any finite 77 . As a 
result, in the extended phase we expect that the asymp¬ 
totic properties can be understood from thermalization, 
i.e., all properties depend only on the energy density of 
the initial state. We compare the properties of our sys¬ 
tem at long times with the thermal state oc 
with P chosen such that ('0o|^|'0o) = Tr(p^j^iL). The 
thermalization arguments are applicable if the subsystem 
of interest (say A) is much smaller than the rest of the 
system {B) as in such situations the subsystem B acts 
as a bath for the subsystem A. Thus we use an unequal 
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FIG. 7. (color online). Asymptotic {t = 10^®) entanglement 
and bipartite fluctuations compared with their values in the 
diagonal ensemble for weak (77 = 1 in a,c) and strong disorder 
(?7 = 10 in b,d). We use an unequal partition {La = 4, Lb = 
8 ). Time averaged entanglement long time after the quench 
agrees with the value predicted from the diagonal ensemble for 
L = 12 in the case of weak disorder (a) while the agreement 
is not good at strong disorder (b). (c-d) The prediction of 
bipartite fluctuations on the other hand is very good for both 
weak and strong disorder. 

partition {La = 4,1/^ = 8) to compare the asymptotic 
entanglement and the thermal entropy. For weak disor¬ 
der (?7 = 1) the system is in an extended phase and we 
find a good agreement between the thermal and asymp¬ 
totic values for both entropy and bipartite fluctuations 
as seen in Fig. (a,c) and their insets. 

For strong disorder we show the asymptotic and ther¬ 
mal values in Fig.[^(b,d) and find a very poor correlation. 
Whereas the thermal values for a given energy density is 
still predicted to be large, both quantities saturate to 
much smaller value in each realization. This is a very 
clear signature of failure of thermalization for localized 
systems. 

D. Comparison to the diagonal ensemble 

MBL systems exhibit an emer^nt integrability in 
terms of local integrals of motion.^^^^^ As a result the 
thermalization picture breaks down, as one needs to take 
into account many conserved quantities, not just the en¬ 
ergy. Th e diag onal ensemble is suitable to handle such 
sitnations) ^^ * ^^ * For a given initial state \'ipo) the diagonal 
ensemble density matrix is defined as 

i 

where |F^i)’s are the energy eigenstates of the system. 
Whereas for a given Hamiltonian the thermal state de¬ 



FIG. 8 . (color online). Gorrelation between entanglement at 
long times after the quench with diagonal ensemble for weak 
(77 = 1 in upper panel) and strong (77 = 10 in lower panel) 
disorder for L = 12. The entanglement is better predicted by 
the diagonal ensemble if the subsystem is much smaller than 
the total system. The dotted lines in the upper panel include 
Page’s correction Sc and is given by Soo = Sd — Sc- 

pends only on the energy of the initial state, the diagonal 
ensemble has more information. By definition the diago¬ 
nal ensemble is obtained by averaging the density matrix 
at all times, hence it can trivially estimate a very long 
time average of any physical observable. However here 
we want to check whether it can predict the properties at 
long times. There is a deficit in entropy when it is mea¬ 
sured foran equal partition if the full system is in a pure 
state.^^^^ To avoid large deficit we use an unequal parti¬ 
tion {La =4,1/^ = 8) while comparing asymptotic and 
diagonal ensemble properties in Fig. We find that the 
diagonal ensemble makes an almost perfect prediction for 
bipartite fluctuations both in the weak and strong disor¬ 
der case, while the prediction for entanglement is good 
only for weak disorder. 

We present the effect of subsystem size in Fig. which 
shows that if the subsystem is much smaller than the to¬ 
tal system the diagonal ensemble makes better predic¬ 
tion for entanglement. For weak disorder the prediction 
is consistent when Page’s correction Sc = dA/{dA S ds) 
is taken into account, dA and ds being the Hilbert space 
dimension of subsystems A and B respectively.!^ 


E. Anderson localization 

So far we have considered the quench properties of a 
disordered system in presence of interactions and the ef¬ 
fects of localization. However a non-interacting system 
(A = 0) localizes in the presence of an arbitrary weak un¬ 
correlated disorder in ID (Anderson localization).!^ The 
localization length will be large at weak disorder and 
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other important difference between the two cases is that 
the entanglement is independent of system size for strong 
disorder in the non-interacting system which is a conse¬ 
quence of the absence of logarithmic growth. 

We confirm that the crossover-like behavior is indeed 
a finite size effect by simulating the quench dynamics in 
much larger systems^ni ^nd observing that the location of 
the local maxima in 6Soo shifts to lower values of r] as 
we increase the system size, see Fig. (d). For a larger 
system size the single particle localization length becomes 
comparable to the system size at a weaker disorder. 

We also note the appearance of different scaling for 
weak and strong disorder from the large system simula¬ 
tion of the non-interacting problem. This different scal¬ 
ing appears to be the limiting behavior with increasing 
system size, i.e., it is independent of system size for large 
enough systems. We speculate that this change in scaling 
is related to the localization length becoming comparable 
to lattice spacing. 


FIG. 9. (color online). Entanglement in the non-interacting 
system, (a) Change in the distribution of asymptotic entan¬ 
glement Soo as a function of disorder ( 77 ). Mean (b) and stan¬ 
dard deviation (c) of Soo versus 77 . (d) Standard deviation of 
Soo vs 77 for large system size using free fermion simulation 
shows different scaling for weak and strong disorder. 


as a result there will be crossover-like behavior as the 
strength of disorder is increased for a finite system. We 
perform global quench simulations of the non-interacting 
system and compare its asymptotic properties to that 
of the interacting system. As noted in earlier studie^^ 
entanglement does not grow logarithmically for any dis¬ 
order in this case but rather saturates after the initial 
rapid increase. The asymptotic distribution shows a bi- 
modal distribution beyond some disorder strength, see 
Fig. (a). For small system sizes, the non-interacting 
system shows a crossover-like behavior and for very weak 
disorder it will effectively be in an extended phase as 
the localization length would be much larger than the 
system size. However even in such situation, there is a 
qualitative difference between the interacting and non¬ 
interacting cases—the fluctuations in saturation value of 
entanglement caused by disorder are very strong for the 
non-interacting case (Fig. as compared to the inter¬ 
acting system (Fig. [^. This is also manifested in the 
standard deviation of entanglement SSoo which decreases 
as disorder goes to zero for the interacting system while 
the opposite behavior is observed for small systems in 
the non-interacting case, compare Figs. (c) and |9](c). 
This difference is seen even for the bipartite fluctuations 
and we expect a similar behavior for all observables. An- 


IV. SUMMARY AND CONCLUSIONS 

To summarize, we have studied the effect of disorder 
on global quench dynamics in a ID spin chain and found 
that the asymptotic behavior of different physical quanti¬ 
ties show signatures of the many-body localization tran¬ 
sition. We first reproduced the result that though the 
mean entanglement after very long time shows a volume 
law for both weak and strong disorder, the value itself de¬ 
creases rapidly as the disorder is increased beyond some 
critical value. More importantly the standard deviation 
of entanglement at large times behaves very similar to 
thermodynamic fluctuations and shows a peak near the 
transition. We find similar behavior for bipartite fluctu¬ 
ations with one important difference, the bipartite fluc¬ 
tuations show an area law in the localized phase and 
its asymptotic values are independent of system size for 
strong disorder. Near the transition the standard devia¬ 
tion of bipartite fluctuations also behave like thermody¬ 
namic fluctuations, a signature that can potentially be 
measured in cold atoms experiments using single atom 
microscopy. We then compared the asymptotic proper¬ 
ties following the quench with canonical ensemble at suit¬ 
able temperature and explicitly showed the breakdown of 
thermalization in the localized phase. This breakdown is 
a result of emergence of quasi-local integrals of motion 
in the system and we make very good predictions of the 
asymptotic bipartite fluctuations once these integrals are 
taken into account by using the diagonal ensemble. We 
finally highlighted the effects of interactions by perform¬ 
ing similar quench study of the non-interacting system. 
We hope that this work would lead to more experimen¬ 
tal efforts to study the MBL phenomena using particle 
number fluctuations. 
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